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1 Introduction: Why SrCu 2 (B0 3 ) 2 ? 

Several properties make SrCu2(B03)2a unique and one of the most interesting quantum magnets [1]. 
This compound is a 2D spin system with a disordered ground state even at very low temperatures and a 
spin gap of about 24 cm -1 (3 meV) in the magnetic excitation spectrum [1, 2]. It has been established 
that the elementary excitations in this system, which consist of a transition from the singlet (S = 0) 
ground state to the lowest energy triplet (S = 1) state, are local, weakly dispersive in the reciprocal space 
while many-'particle' magnetic states are more mobile, i.e. show more dispersion. The strengths of the 
relevant magnetic interactions place this compound close to a quantum critical point (QCP) separating the 
gapped phase from a gapless state having long range magnetic AF order. Moreover data in high magnetic 
fields show plateaus at commensurate (1/8, 1/4 and 1/3) values of the saturation magnetization [2, 3, 
4]. The plateau states can be thought of as crystalline arrangements of magnetic moments separating 
regions of continuous rise in magnetization, these latter regions allowing for an interpretation in terms 
of Bose-Einstein condensation of triplet excitations [5]. It has also been suggested [6] that doping in 
this system (regarded as a Mott-Hubbard insulator) may lead to a superconducting phase mediated by 
antiferromagnetic (AF) fluctuations, a mechanism similar in spirit to one of the scenarios proposed for 
the high T c cuprates [7]. 

The first synthesis of SrCu2(B03)2 was achieved in 1991 by Smith and collaborators [8] but the authors 
did not elaborate on the magnetic properties. This compound was rediscovered in 1999 by Kageyama et 
al. [2], who also pointed out the outstanding magnetic properties and the importance of SrCu2(B03)2 in 
the physics of low dimensional quantum magnets. The crystal belongs to the tetragonal symmetry and it 
has a layered structure in which Cu(BOs) units are separated by closed shell Sr 2+ atoms, see Fig. 1. At 
T = 395 K the system undergoes a 2 nd order phase transition from the space group 14/ mem to J42m on 
cooling from high temperature. In the IA/mcm phase the planes containing the Cu atoms are flat and 
they form mirror symmetry elements. Below 395 K the transition can be intuitively understood as the 
buckling of the Cu planes which lose their mirror symmetry. 

The magnetic properties of SrCu2 (603)2 are determined by the S = 1/2 spins sitting on the Cu 2+ 
sites. The crystal structure imposes a distribution of the magnetic moments in weakly coupled 2D layers 
which define the (ab) planes. In each of these sheets the spins form orthogonal dimcr lattices and the 
c-axis is perpendicular to them. The in-plane magnetic interactions can be described in an effective spin 
model by taking into account the magnetic exchange between nearest neighbor Cu spins bridged O atoms 
forming a 147° Cu-O-Cu bond which points toward an AF coupling denoted by Ji, see Fig 2. Besides 
this interaction one can consider the next nearest neighbor super-exchange taking place via the BO3 
complexes. Curie- Weiss type fits of the high temperature magnetization data suggest that this coupling, 
denoted by J2, is also AF. Indeed, as we will discuss later, these two terms seem to capture many aspects 
regarding the magnetic properties of SrCu 2 (60-3)2 ■ The addition of a weak inter layer exchange J' to 
the above two terms constitute the starting point for treating the 3D spin dynamics in this compound. 

The schematic of the magnetic lattice is shown in Fig. 2 where the nearest neighbor AF bonds are 
represented by a solid line while the next nearest neighbor interaction is shown by dashed lines. Note an 
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Fig. 1. 3D view of the SrCu2(BC>3) 2 crystal showing the layered structure along the c-axis. Above 395 K the 
Cu(BOs) planes are flat, mirror symmetry elements. 

important property of this lattice: ignoring the solid bonds (the J\ interactions), the topology of the spin 
structure, determined by the J2 bonds, is equivalent to that of a 2D square lattice. The corresponding 
Hamiltomian for one Cu(B03) plane reads: 
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Here i and j are nearest neighbor (NN) or next nearest neighbor (NNN) Cu sites. 

Why is a Raman study of this compound interesting? Several magnetic properties, features of the 
excitation spectrum or the occurrences of the magnetization plateaus, can be explained by taking into 
account the J\ and Ji interactions. As will be discussed in the review of the main theoretical results, 
this compound is believed to be characterized by a ratio x = J2/J1 close to 0.7, which is the value 
around which the system ground state changes. One reason to study this compound is that an accurate 
determination of this ratio is still missing and this is important due to the proximity to the QCP. In 
this regard the recent observation of a mode below the spin gap brings into question the quantitative 
determination of the AF exchange parameters of the system since the existence of such excitation has 
not been predicted by theory. Another reason is brought about by a different set of questions related to 
the way the external radiation field couples to magnetic excitations. Magnetic modes in the triplet sector 
have been probed spectroscopically, besides inelastic neutron scattering (INS) [2, 16], also by electron 
spin resonance (ESR) [17, 18] and infra-red (IR) [19, 21] absorption experiments. The nature of the 
mixing interactions originating in spin-orbit coupling which usually allows transitions from the singlet 
(S = 0) to excited (S = 1) states is still to be understood. Our approach, which is a study of collective 
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Fig. 2. Left: The (06) planes and crystal axes notation. Right: The effective magnetic lattice with the unit cell 
of the dimer structure being shown by the short dashed square. The circles are the 5 = 1/2 Cu spins, the solid 
and long dashed lines represent the intra-dimer (Ji) and inter-dimer (J2) AF superexchange interactions. Below 
395 K [9] the vertical and horizontal dimers become slightly non-coplanar, see Fig. ??. 
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magnetic excitations in terms of symmetry, resonance and coupling mechanisms in external magnetic 
fields, is illuminating in this respect. 

This chapter will be focussed on the low temperature properties of SrCu 2 (603)2 ■ The discussion of 
the low temperature phononic excitations as seen in Raman spectra in the next part will be followed by a 
review of experimentally found magnetic properties with an emphasis on spectroscopic techniques. Then 
we will present a description of the basic properties emerging from the Hamiltonian (1) which describes 
a Shastry-Sutherland lattice and the structure of the magnetic excitation spectrum. The effects of other 
interactions, in particular of the antisymmetric Dzyaloshinskii-Moriya (DM) terms, will be described. 
Then we will present our low temperature Raman data and discuss our results in reference to the open 
questions mentioned above. 

2 Low Temperature Phononic Spectra in SrCu 2 (B0 3 )2 

The analysis of the low temperature lattice dynamics is important for several reasons. In general this is 
because there are many examples of low dimensional crystalline compounds which undergo transitions 
to phases in which the appearance of a spin gap in the magnetic excitation spectrum is accompanied by 
real space lattice symmetry breaking due to spin-phonon coupling. The lowering of the crystal symmetry 
may involve newly allowed phononic modes which could be checked directly in the Raman spectra. 
In the particular case of SrCu2 (603)2 spin-lattice interaction has been suggested to be relevant to the 
magnetic dynamics at low temperatures and/or high magnetic fields. In order to explain the selection 
rules of the transitions seen in IR absorption, spin-phonon induced antisymmetric DM spin interactions 
have been invoked [21]. It has been argued that this coupling will induce virtual phonon transitions 
which will instantaneously lower the crystal symmetry, allowing for non fully-symmetric effective singlet 
triplet mixing terms. On the other hand the spin-lattice coupling was taken into account in order to be 
able to describe the spin density profile at high fields in the magnetization plateaus states. NMR data 
brought evidence for the broken translational symmetry at the 1/8 plateau by the existence of at least 11 
nonequivalent Cu sites [4] . The role of the phonons coupled adiabatically to the spin degrees of freedom 
in this case is to lift the degeneracy of the ground state, picking a state with a certain magnetization 
texture and allowing the ground state magnetization to be a defined quantity. 

So far the study of phononic excitations has been especially focussed in relation to the structural 
phase transition at T = 395 K [9, 10]. The lattice soft mode of the 2 nd order transition from the 1 4/ 'mem 
to the /42m group can be seen in Fig 3. The phonon which condenses belongs to the Bi u representation 
of the higher symmetry group and, in terms of Cu atoms, it involves mainly an alternate displacement 
along the c-axis of the nearest neighbor dimers, see Fig 1. 

A symmetry analysis of the phononic excitations in SrCu2 (603)2 is done in the following. The unit cell 
contains 4 formula units and a total of 44 atoms. The number of k = modes is given by 3 x 44/2 = 66, the 
factor of 2 coming from the fact that the unit cell is body centered. In this case, the symmetry lowering at 
395 K does not change the number of atoms in the unit cell. The associated point group of the I A/ mem 
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Fig. 3. The fully symmetric soft mode below the IA/mcm - /42m transition from Ref. [10] seen in cc polarization. 
Note the strong two-phonon excitation around 120 cm - . 
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Table 1. Symmetry of the phononic excitations of SrCu2(BC>3)2 in the high (74/mcm, point group D±h) and 
low (742m, point group Did) temperature phases. The first column contains the atom types and the rest of the 
columns correspond to irreducible representations of the point groups. 
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space group is D±h while 7?2d corresponds to 742m. The analysis is based on the tables in Ref. [23] and 
uses the site symmetry approach. Oxygen atoms mediating the intra-dimer superexchange (01) and the 
Oxygens allowing for the inter-dimcr superexchange (02) occupy a different symmetry positions. In the 
high temperature phase the Cu, 01 and B atoms have C' 2vl 02 atoms have C' s and Sr atoms occupy D4 
site symmetries respectively. Table 1 summarizes the number of modes corresponding to each atom both 
above and below the transition. The part of the table related to the 742m group can be easily inferred 
from the analysis of the high temperature phase. One has to drop the 'u' and 'g' indices (corresponding to 
odd and even modes) because the inversion symmetry is lost and by using compatibility tables which show 
that the representations Ai u , A 2u , B\ Ul B 2u of the D^h point group become By, B%,Ax, A 2 representations 
(in this order) in D 2 d- Both the E u and E g representations remain double degenerate. 

Inversion symmetry breaking results in the fact that certain phonons, which were dipole active, become 
Raman allowed. The appearance of several new modes below 395 K is shown in Fig. 4 where data from 
Refs. [9, 10] is reproduced. Notable are the shoulders appearing around 160 and 220 cm" 1 below 395 K 
in ca polarization (which probes double degenerate modes with E symmetries) suggesting either an 
almost degeneracy of E u and E g modes in the high temperature phase or, even more interestingly, the 
possibility that the symmetry of the low temperature phase is lower than what has been inferred so far 
and accordingly, the modes corresponding to higher dimensional representations become non-degenerate. 
[9, 10] is reproduced. 

We show in Fig. 5 our low temperature Raman data taken with in plane polarizations and using 
u>l = 1.92 eV laser excitation energy. In the D 2 d group the (RR), (RL), (aa), (ab), (a' a') and (a'b') 
probe A\ + A 2 , B\ + B 2 , A\ + B\, A 2 + B 2 , Ay + B 2 and A 2 + B\ symmetries respectively. The modes 
below 60 cm" 1 are not indexed since they are magnetic and will be discussed in the next section. The 
modes are sharp and they are sitting at this low temperature on a flat background with almost vanishing 




Fig. 4. Phononic data from Ref. [9] (left) and Ref. [10] (right) showing in plane and out of plane polarized spectra 
above and below the structural transition. Newly allowed Raman phonons are marked by arrows. 
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Fig. 5. Raman data at T ~ 3 K in six polarizations taken with the 1.92 eV laser excitation energy. The energy 
and symmetry of the modes above 60 cm -1 are shown in the figure. 



intensity. The 60 cm -1 mode corresponds to the soft mode of the structural transition. At 121.8 cm -1 
we observe in the Ai channel the two phonon excitation seen also in (cc) polarization, Fig. 3, and very 
close to it another sharp mode with B\ symmetry. We observe, similarly to the spectra shown in Fig. 4b, 
several pairs of modes having very similar energies. For example doublet structures are observed around 
284 cm -1 where we see a pair of Ai and Bi excitations and also two modes having Bi and B2 symmetries 
are found around 320 cm -1 . 

One way to explain this behavior is to follow up the suggestion in Ref. [10] and assume that in the 
high temperature phase there are phonons which are odd and even with respect to inversion but very 
close in energy and to try to identify them by looking at similar atomic vibrations corresponding to 'u' 
and 'g' representations respectively. The (Ai,Bi) group around 284 cm -1 in Fig. 5 would correspond 
in this scenario either to a group of (Ai ff ,Ai„) or (Bi 9 ,Bi„) in the high temperature phase since Ai u 
representation becomes Bi and Bi u representation becomes Ai at low temperatures. Similar reasoning 
would suggest for the (Bi,B 2 ) group around 320 cm -1 that it originates either from a pair of (Bi ff ,A 2u ) 
or (B2 S ,Ai„) modes above 395 K. We performed a symmetry analysis of the k = atomic vibrations. 
The conclusion is that this approach does not provide an a priori reason for the quasi-degeneracy and 
one has to perform a quantitative normal mode energy calculation by using appropriate inter-atomic 
clastic constants. There is a simple way to see why this is true. Looking at the character table of the D±h 
group it can be noticed that the even modes are symmetric with respect to the mirror symmetry in the 
Cu(B03) planes while the odd modes are antisymmetric with respect to this symmetry operation. This 
means that the 'u' phonons in the high temperature phase correspond to vibrations of the atoms along 
the c-axis while the 'g' modes consist of in-plane movements. Due to the different oscillation pattern 
it is hard to explain the closeness of phononic energies at this qualitative level. We note that one can 
easily find vibrations corresponding to different group representations which involve similar oscillations 
at the 'molecular' level (for instance groups of atoms forming the Cu-0 spin dimer structure or pairs 02 
atoms between nearest neighbor dimers) and having a certain inter-molecular phase pattern, but they 
do not correspond to the experimentally observed symmetries. We suggest that, remaining within the 
conclusions of X-ray studies which have not found evidence for additional crystallographic changes at low 
temperatures, good candidates for understanding this intriguing behavior are provided by the inter-dimer 
BO3 molecular complexes whose rotations as a whole around the a' and c-axes or whose in and out of 
the plane translations may turn out to be similar in energies. 
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Fig. 6. Temperature dependence of the magnetization (a) and inelastic neutron scattering results as a function of 
in plane wavevectors at low temperatures in SrCu2(BOs)2 (b) from Ref. [2]. The solid line in (a) shows the result 
of a fit from numerical calculation using x = J2/J1 = 0.635, Ji = 59 cm -1 (7.32 meV) and also an interlayer 
coupling J3 = 5.5 cm -1 . The interactions Ji and J2 alone provide a good description of the low temperature data, 
including specific heat, but a finite J3 was necessary to explain the high temperature behavior of the magnetization 
where a Weiss temperature 9 = -92.5 K was obtained from the fit between 160 and 400 K [2]. 



3 Magnetic Properties of SrCu 2 (B03)2 
3.1 Experimental and Theoretical Reviews 

Experimental review — In this part we discuss data which relate to the most interesting properties 
and set the relevant energy scales of SrCu2(B03)2 ■ Magnetization and INS data provided for the first 
time evidence for the existence of a gapped phase in this compound [2]. Fig. 6a shows that there is a 
drop in x below about 15 K suggesting a gapped phase. Another important thing pointed in the inset of 
this figure is the strong suppression of the magnetization peak around 20 K compared to the prediction 
of a simple dimer model which points towards the importance of other magnetic interactions. Neutron 
scattering, Fig. 6, probed the excitations out of this phase and confirmed the existence of a gap of about 
3 meV (24 cm -1 ) and found additional excitations around 5 and 9 meV. Notable is the flat dispersion 
(less than 6%) of the gap branch seen at 3 meV as a function of in-plane wavevectors meaning that 
the lowest excitations are very local. The dispersion of the higher energy branches is more pronounced 
suggesting more mobile excitations in the multi-triplet channels. The inability of the dimer model to 
describe the experimental data suggests that inter-dimer interactions are important and points towards 
a very interesting physics of frustration in this system. 

Besides INS and magnetization, ESR and IR data in magnetic fields confirmed the magnetic nature of 
these excitations. These experiments, along with high resolution INS data display a rich internal structure 
of the magnetic excitations and various selection rules for in and out of the plane applied magnetic fields, 
see Fig. 7. In particular, the gap multiplet is shown to be made of 6 branches, which is expected since the 
unit cell contains 4 spins and therefore exciting a triplet on each of them while the other one is in a singlet 
state will give a total of 6 excitations. Out of these 6 modes, 4 are seen to be symmetrically disposed 
around the gap value A = 24.2 cm -1 . The observed splitting suggests that besides the superexchange 
interactions there are other low energy scale interactions which have to be taken in to account. 

Theoretical review — This section discuses the salient properties of the magnetic structure of 
SrCu2(B03)2 and the excitations of the Hamiltonian (1). Among them, the ground state properties, the 
localization of the elementary triplets leading to the dispersionless magnon branches from Fig. 6b, the 
existence of strongly bound two-triplet states as well as the presence of additional spin orbit couplings 
generating the fine structure seen in Fig. 7. 
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In the approximation given by Eq. (1) the wavefunction given by the direct product of singlet dimers 
is always an eigenstate of the system and it can also be shown that it is the ground state for a continuous 
set of parameters x = Jij J\. For the proof of the eigenstate one has to take into account only the second 
term in (1) since the dimer product is obviously an eigenstate of the AF coupled independent dimers. 
The application of the J2(Si • S3 + S2 • S3) on the \s\i > ®|s34 > vanishes because the operator (Si + S2) 
applied to the singlet state |si2 > vanishes. This is essentially due to the different parity of the singlet 
and triplet states with respect to the 1 <-> 2 exchange and the fact that on the orthogonal dimer lattice 
the Hamiltonian conserves the parity. This remains true even if the interlayer coupling is considered [1] . 

For x < 0.5 one can show that the singlet dimer product, \tp >, is indeed the ground state. The 
Hamiltonian (1) can be written as H = J^f* hi where Nt is the number of triangles of the type formed 
by the spins 1, 2 and 3 in Fig. 9 and h, = (Ji/2)(Si • S 2 ) + J 2 (Si • S 3 + S 2 • S3). The ground state of 
each hi is e* = — 3/8Ji if x < 0.5. Accordingly, denoting the true ground state of H by \<j> >, one has 
Eg = < 4>\H\4> > — X^i ' < 4>\hi\4> > > — 3/8JiA/j because \<f> > is a variational function for hi. Considered 
as a variational wavefunction for H and taking into account that the action of the second term in Eq. (1) 
on \ip > is identically zero, one obtains immediately that < ip\H\ip >= —3/4 J\Nd = —3/8 J\Nt because 
the number of dimers, Nd, is half of that of the triangles. On account of the variational principle, the true 
ground state energy satisfies E g < —3/8J\Nt. From the two inequalities one obtains that E g = —3/8J\Nt 
and so the product of singlet dimers is indeed the ground state. The underlying reason for these beautiful 
properties is the fact that the magnetic lattice of SrCu2 (1303)2 is a realization of a 2D Shastry-Sutherland 
model which up to now has been discussed only at a theoretical level. 

Turning to the question of the flat dispersion of the excitations in the one-triplet sector one can 
consider in Fig. 9 the case of two nearest neighbor orthogonal dimers having an excited triplet state on 
the vertical bond and a singlet on the horizontal one. The propagation of the triplet on the horizontal 
bond by the NNN term in Eq. (1) is only possible if a triplet is left behind on the vertical dimer. This is 
due to the reflection symmetry of singlet and triplet excitations with respect to a mirror plane parallel to 
the horizontal bond, Fig. 9a and the fact that the Hamiltonian of the system must be a fully symmetric 
operator. If a triplet is left behind, then the hopping of one triplet to a neighboring dimer is very restricted 
and is possible only by virtually forming closed paths of triplets, the smallest of these paths involving 
three adjacent dimers. As a result, one triplet hopping appears only in the 6 th order of perturbation 
theory showing that these excitations are very localized in real space and explaining the flat k dispersion 
seen in Fig. 6. The motion of two triplets is different however. It has been shown [11] that in this case 
correlated hopping processes can occur and it was found that two-particle hopping appears in 2 nd order 
perturbation theory. This explains the larger dispersion seen by INS for the excitations around 5 meV. 

We saw that the products of dimer singlets is the ground state of H at least for x < 0.5. What is the 
value of x above which this is no longer true? That there should be a finite value is obvious from the fact 
that in the limit x — > 00 the lattice is topologically equivalent to a 2D square lattice which has long range 
order at T — K and a spin- wave approximation is more appropriate. Theoretical work [1, 11, 12, 14] 




H (T) H (T) 

Fig. 7. Neutron (Ref. [16]) and ESR (Ref. [17]) data at low temperatures as a function of magnetic fields parallel 
(a) and perpendicular (b) to the c-axis. 
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Fig. 8. Predictions for the two-triplet bound states in the singlet (left) and triplet (right) of a perturbative 
analysis of the Shastry-Sutherland Hamiltonian of Eq. (1) from Ref. [11]. A, i~2, A, A and Zg correspond to j4i, 
Si, B2, ^2 and Z? in the notation of Table 1. 

shows that below a ratio x — Jij J\ « 0.7 the ground state remains the same and the system has a finite 
spin gap, A to the lowest excited S = 1 state. The spin gap is equal to Ji for a; = but with increasing 
this ratio A gets renormalized down due to many body effects, see Fig. 8. At high values of x the system 
has long range order, other possible intervening states separated by QCP's have been proposed to exist 
around x = 0.7. The gap renormalization as a function of x is shown in Fig. 8. As shown here and also 
in the parameters used to fit magnetization data in Fig. 6, SrCu2(B03)2 is closed to a QCP having a gap 
renormalized significantly, to a bit more than 50% of the 'bare' value given by J\ sa 85 K (59 cm -1 or 
7.3 meV), see Ref. [1]. 

Another characteristic of this magnetic lattice seen in Fig. 8 is the existence of bound states in 
the two-triple sectors, see Fig. 8. These are states which have an energy below the inset of the two- 
magnon continuum starting at 2A. Many of these have been predicted along with their symmetries. 
The vertical line in Fig. 8 showing the position of SrCu2 (603)2 in the phase diagram was inferred from 
the experimentally found values of the spin gap A = 24 cm -1 (3 meV) and the observation by Raman 
scattering of a strong and sharp resonance (attributed to collective 5 = two-triplet bound state) 
around 30 cm -1 [20]. We will discuss such excitations in the next section devoted to the analysis of low 
temperature Raman data in SrCu2(B03)2 ■ 

Before that, another observation in connection to the experimental results shown in Fig. 7: there it is 
seen that the 6 branches of the gap multiplet around 24 cm -1 are split in three pairs of doublets. Cepas 
and collaborators proposed that this is due to the existence of inter-dimer DM interactions which have 
a direction parallel to the c-axis, see Fig. 9. The proposed Hamiltonian (for the spins in the unit cell) 
reads: 

h c = d c 13 (S 1 x S 3 ) + d c 14 (Sr x S 4 ) + d c 23 (S 2 x S 3 ) + d c 24 (S 2 x S 4 ) (2) 

The DM vectors satisfy d c 13 = d c 24 = — d c 14 = — d c 23 due to the crystal symmetry, in particular due 
to the existence of mirror planes orthogonal to the dimers as shown in Fig. 9 and the existence of C2 
rotation axes parallel to the c-axis and passing through the middle of the dimers. It was found that 
this interaction reproduces the behavior in magnetic fields perpendicular and parallel to the (ab) plane, 
which is shown by solid lines in Fig. 7. The energy of the upper and lower pairs is given by A ± d c 
and they remain degenerate, see also Fig 11c where we show symmetry analysis results. Importantly, the 
inter-dimer DM interactions parallel to the c-axis are allowed both above and below the structural phase 
transition at 395 K. Other DM terms, for instance the red arrows in Fig. 9, the are not allowed above T c 
because the Cu(B03) plane is a mirror symmetry element. 
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Fig. 9. Left: Cartoon with two nearest neighbor dimers discussed in the text related to the exact ground state 
and hopping of triplet excitations, m denotes a mirror plane; the singlet and triplet states on the 1 — 2 dimer will 
have different symmetries with respect to this reflection operator. Right: The unit cell of the magnetic lattice from 
Fig. 2 with the inter-dimer Dzyaloshinskii-Moriya term suggested in Ref [16] which is parallel to the c-axis. The 
arrows perpendicular to each dimer correspond to our proposed antisymmetric intra-dimer interaction leading to 
singlet-triplet mixing. 



3.2 Magnetic Raman Scattering Results in SrCu2(BOs)2 

Here are some experimental details related to our data discussed in the following. The spectra were taken 
from the (ab) single crystal surface in a backscattering geometry. We used an incident power of about 
0.6 mW focussed on a 100 /i diameter spot. The crystallographic axes orientation was determined by X- 
ray diffraction. The data in magnetic fields, taken at a sample temperature of about 3 K, were acquired 
having the continuous flow cryostat inserted in the horizontal bore of a superconducting magnet. We 
used the u>l = 1-92 and 2.6 eV excitation energies of a Kr + laser and a triple-grating spectrometer for 
the analysis of the scattered light. The data were corrected for the spectral response of the spectrometer 
and detector. Polarization configurations are denoted by (ei„e out ) where these two vectors are along the 
polarization direction of the incoming and outgoing photons. Circular polarizations are denoted by (RR) 
and (RL) where ei„ iOUt = (d ± ib)/\/2. 

In Fig. 10 we show six low temperature Raman spectra in zero applied field and using a laser frequency 
lul = 1-92 eV. The symmetries probed by each polarization in the tetragonal group are shown in the 
legend. The mode seen at 60 cm -1 is the soft mode of the structural transition, also seen in Fig. 3, which 
belongs as expected to the fully symmetric representation. Three strong features are seen in the spectra at 
23, 29 and 38 cm -1 and they are denoted by Tib, Si and T2. These excitations belong to the A2, Bi and 
A2 representations respectively. Besides these three modes, we observe several other weaker excitations. 
Among them we see a set of three A2 symmetric modes denoted by T3, T4 and T5. One can also observe 
the presence of other very weak feature at 25.5 cm -1 in all polarization configurations except (RR). This 
suggest that at this frequency there are two quasi-degenerate excitations and that they belong to the 
Bi and B2 representations which would justify their observations in five out of six scattering geometries. 
The symmetry analysis discussed later confirms indeed the above assumption. The modes are denoted 
by Ti e and Tif ■ A summary of the zero field energies and their experimentally observed symmetries can 
be found in Table 2. The energy at which we observe these modes, the comparison with the INS, ESR 
and IR data, part of which are shown in Figs. 6 and 7 as well as our data in magnetic fields confirm the 
magnetic origin of these modes and the predominant S = 1 character of the 'T' modes. The Tib, Ti e 
and Tif excitations modes seem thus to belong to the spin gap multiplet while T2, T3, T4 and T5 would 
correspond to multi-particle S = 1 channels. 

The zone center elementary excitations are generated by the spins within the magnetic unit cell. If 
the picture of real space localized elementary triplets is true, then one expects that the analysis of the 
4-spin cluster forming the unit cell is able to predict correctly the experimentally observed symmetries 



10 



A. Gozar and G. Blumberg 



o 



co 
E 

CO 
DC 



(RR):A i+ A 2 / 

-^(ab):A 2 + B 2 

^(a'b'):A 2 +B i 
. — (aa) : A, + B 1 


IT 

1b 

c 

a 1 

1 


T 

2 


- 


_^(RL):B i + B 2 I 


T T 

1e 


1 




—a— (a'a') : A + B f 


T 

n J 
V J 


A T T T 

/\ 3 4 5 


- 

K. 





























10 20 30 40 50 60 



Raman shift (cm" 1 ) 

Fig. 10. Zero field Raman data in SrCu2(B03)2 taken with u>l = 1-92 eV excitation energy at T — 3 K in six 
polarizations. The legend shows the tetragonal symmetries probed in each scattering geometry. 



of these excitations. Accordingly we calculated the energies and symmetries of the excitations generated 
by the magnetic cluster shown in Fig. 11a and the Hamiltonian of Eq. (1). The total number of states, 
16, consist of two singlet, 5 = 0, states, three triplet, 5 = 1, states and one quintuplet, 5 — 2 state. The 
energies are shown in Fig. 11a as a function of x — Ji/ J\ and one can observe their linear variation with 
the this ratio. The points at x = and x = 1 can be easily understood. The first one corresponds to 
two independent AF coupled by J\ on the bonds 1 <-» 2 and 3 <-> 4 so the energies will be just the sum 
of those corresponding to singlet and/or triplet states sitting on each dimer separately. For the second 
point one can see that [H, S tot ] = so the energies will be, up to a constant factor, equal to the total 
spin eigenvalues 5(5 + 1) (in units of h 2 ). Starting with four 5=1/2 spins one can have for the total 
spin the values 5 = 0,1 and 2. 

The cluster in Fig. 11a belongs to the Did group which is, as it should be, the point group associated 
with the space symmetry of the crystal. We performed the group symmetry analysis of the eigenstates 
by decomposing the representation obtained starting with the canonical spin basis (the direct product of 
I j> and I i> spinors representing each Cu atom) into the irreducible representations of the Did point 
group. Each symmetry element of the Did group involved the action on spinor state in the origin of the 
coordinate system as well as the corresponding permutation of spin indices. Using projection operators 
we block diagonalized the Hamiltonian and numerical analysis (including the DM terms shown in Fig. 9 
and discussed in more detail later) allowed us to identify the energy and the predominant spin character 
of each eigenstate. 

Indeed, as it is shown in the right-most column of the Table 2 and Fig. 11c, the symmetries of the 
observed one-triplet excitations correspond to the results of group theory analysis. The latter predicts 
for the six elementary triplets the Ai, A2, Bi, B2 and E symmetries, see Fig. 11. The E modes, Ti c 
and Ti,j, cannot be observed in zero field when the light propagates parallel to the c-axis because they 
are accessible only in (ca) or (cb) polarizations. The fully symmetric Ti a mode which, within the spin 
model including J\, Ji and the out-of-plane inter-dimer DM interaction proposed in Ref. [16], should 
be degenerate with the strong Ti& (A2 symmetric) excitation at 22.8 cm -1 is also not observed. This is 
most probably due to a much weaker coupling to light in this symmetry channel. The analysis shown in 
Fig. 11 shows that the pair of Ti e and Ti / modes at 25.6 cm -1 are to be seen in the Bi and B2 symmetry 
channels which explains the weak feature seen in Fig. 10 in all polarizations except the (RR) scattering 
geometry. The observation of the Ti e (Ti/) modes with Bi (B2) symmetries 2.8 cm -1 above the A2 
symmetric mode T± e allows the determination of the magnitude of the inter-dimer interaction d c (see 
Fig. 11) and also of its absolute sign because a sign change will interchange the position of these modes 
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(a) (b) (C) 




Fig. 11. (a) The cluster used for the symmetry analysis of the k = modes generated by the four spins in the 
unit cell of Fig. 9. (b) The eigenvalues of the Hamiltonian corresponding to the 4-spin cluster in (a) taking into 
account the Ji and J2 terms as a function of the ratio x = J2/J1. (c) Results of the symmetry analysis for this 
cluster in the Z?2d group. The energy levels correspond to a value x = 0.3. On the right we show the energy 
splittings and the absolute and relative to the ground state (in parenthesis) symmetries of the 16 magnetic modes 
when a finite inter-dimer DM interaction d c is present. A finite intra-dimer DM term d a b will further split the 
modes which belong to the one-dimensional representations, see Fig. 9. 

around the gap A. Finite intra-dimer Dzyaloshinskii-Moriya interactions (red arrows in Fig. 9) will split 
slightly the energies of the states in Fig. 11 which belong to the one dimensional representations. These 
effects will be discussed in more detail later where we analyze quantitatively the behavior of the spin 
excitations in magnetic fields applied parallel and perpendicular to the c-axis. 

We turn now to the discussion of two-triplet states. Besides the singlet ground state and six one-triplet 
states, the symmetry analysis of the 4-spin cluster shown in Fig. 11 predicts the following: one S = 
two-triplet bound state to be observed in the Bi g channel, three branches with Ai and double degenerate 
E symmetries which belong to another bound S = 1 excitation and five branches of a quintuplet (S = 2) 
state to be accessed in the Ai, A2, Bi and E channels. Their symmetries, along with the energies for the 
particular value x = 0.3 chosen as an example, are shown in Fig. 11c. We observe that due to symmetry 

Table 2. Collective spin excitations in zero field: notation, the predominant spin character, Stot and 2 projection 
S z , energy and transition symmetries as observed experimentally and predicted from the 4-spin cluster in Fig. 1 
corresponding to k = excitations. 



Mode 


Stot (S 2 ) 


Energy 


Symmetry 








Experiment 


Group Theory 


Tla 


1 (±1) 


22.8 




Alg 


Tu 


1 (±1) 


22.8 




A 29 


Tic 


1 (0) 


24.2 




E 


T w 


1 (0) 


24.2 




E 


T le 


1 (±1) 


25.6 


Big 


Big 


TV 


1 (±1) 


25.6 


B 2o 


B 2 g 


Si 


0(0) 


28.9 


Bi 9 




T 2 


1 


37.5 


A 29 




T 3 


1 


40.8 


A 2g 




S 2 





43.0 


Big 




T 4 


1 


44.5 


A 29 




S3 





45.2 


Big 




T 5 


1 


50.9 


A 29 
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reasons none of the A 2 symmetric modes from T 2 to T 5 having energies higher than 30 cm -1 qualify for an 
interpretation as triplet bound states generated within the 4-spin cluster. This is consistent with the fact 
that larger cluster sizes are necessary in order to capture the more delocalized nature of these excitations 
which implies that they have contributions from the different parts of the Brillouin zone. The fact that 
the existence of the strong A 2 symmetric bound triplet state at an energy 1.55 -A = 37.5 cm -1 has 
not been predicted by high order perturbative analysis [11] suggest that other spin interactions have to 
be taken into account in order to explain the excitation spectrum. Apparently symmetry considerations 
would allow the 28.9 cm -1 feature denoted by Si in Fig. 10 to be interpreted as the singlet bound state of 
two triplets within a unit cell. As we show in the following section, this 28.9 cm -1 mode does not shift in 
external fields, which is compatible with a collective singlet excitation as discussed in [20], but suggests 
that its internal structure is not the one derived from the 4-spin cluster. 

In Fig. 12, using the same mode notations, we show the influence of an external magnetic field applied 
parallel and perpendicular to the c-axis on the low temperature Raman spectra from Fig. 10. The relevant 
aspects are the following. In panel (a) we observe the splitting of the Ti a and Tit, modes in magnetic 
fields B || c, the B = 1 T showing that the A 2 mode (T^) present in zero field disperses upwards with 
increasing the magnitude of the field. Dashed lines in this figure mark the dispersion of the much weaker 
modes T3, T4 and T5. In Fig. 12b one of the E modes becomes Raman active due to symmetry lowering 
for B _L c configuration and we observe three dispersing branches of the gap multiplet. Fig. 12c shows 
that the Bi symmetric excitation at 28.9 cm -1 does not change its energy with field, only a very small 
negative shift of the order of 0.5 cm -1 from to 6 T is seen because of the crossing with the upward 
dispersing gap branches seen in (RR) polarization. Panel (d), which is a zoomed in region of Fig. 12b, 
shows that several modes become Raman active in finite fields B _L c around 38 cm -1 where the T 2 
excitation lies. The internal structure of this higher energy multiplet is composed of modes dispersing 
up, down or independent of magnetic field. Remarkable is their similar selection rules and dynamics in 
magnetic fields of the collective modes around 38 cm -1 and 24 cm -1 . The emergence in finite fields of 
several strong modes in the spin gap region precludes the observation of the weak Ti e and Ti / modes 
from Fig. 10. 

Fig. 13 summarizes the magnetic field dependencies of the energies and spectral weights of the most 
intense Raman excitations. The symbols in Fig. 13a-d correspond to experimental data, dashed lines 
are guides for the eye while the solid lines are results of a numerical diagonalization of a 4-spin cluster 
using the same set of parameters. Taking into account that the 4-spin cluster neglects many-body gap 
rcnormalization effects (see Fig. 8) leading to a singlet-triplet energy independent of x = J 2 /Ji as well 
as the fact that when using periodic boundary conditions there is an effective doubling of the J 2 and 
inter-dimer DM interactions from Eq. (2), we chose the following values: J\ = A = 24.2 cm -1 which is 
the value of the spin gap, see the Table 2; x — 0.556 from the ratio of the energies of the sub-gap mode 
at 21.5 cm -1 [18, 21] with respect to the gap A (see Fig. 11); an inter-dimer DM term parallel to the 
c-axis, d c — 1.4 cm -1 , which produces the splitting of the Ti a; f, and Ti e j branches from 24.2 cm -1 (our 
value is consistent to the one proposed in the literature [16]); finally, from the magnetic field value around 
which the intensities crossing in Fig. 13 takes place, the value chosen for the intra-dimer interaction was 
d a b = 2.66 cm -1 . 

The term containing the intra-dimer DM interaction is proposed by us in order to explain the coupling 
and selection rules with magnetic fields. Symmetry considerations impose in the I42m group for the 
directions of the in plane DM vectors the ones depicted in Fig. 9 with the Hamiltonian reading: 

h ab = di 2 (S! x S 2 ) + d a f(S 3 x S 4 ) (3) 

The interaction terms responsible for the coupling to the external photon field and which were used for 
the calculation of intensities in Fig. 13c-d are discussed in the next part. We also plotted (filled squares 
and dashed lines) the field dependence of other higher energy modes observed in Fig. 12. 

We remark an overall qualitative agreement for this choice of parameters which suggests that the intra- 
dimer interaction has to be taken into account in the spin Hamiltonian. The agreement is quantitative as 
regards the energies and the intensity variations for the B 1 c case. The term described in Eq. (3) plays 
a crucial role in obtaining a finite coupling to the excited S = 1 triplets which is not realized by the d c 
terms in Eq. (2). Regarding the apparent degeneracy of the Ti a and Tu, as well as that of the Ti e and 
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Fig. 12. Magnetic field dependences of the magnetic excitations at T = 3 K using the ujl = 1.92 eV excitation 
energy in the following geometries: (a) (RR) B || c. (b) (RR) B _L c. In (c) the (RL) polarized data is shown for 
and 6 T magnetic fields for both B || c and B _L c. In (a) and (b) the vertical shift is proportional to the magnetic 
field difference with respect to the 9 T spectrum and the dashed lines are guides for the eye. 

Tif modes although they belong to different representations, we note that these two groups of excitations 
are split by h a b but the splitting is very small, of the order of d^ h j A 0.25 cm -1 . The largest discrepancy 
between the experimental data and the calculation is seen in Fig. 13. One aspect in this regard is that 
the value d a b had to be chosen greater than that of d c . This is intriguing because the d c term is allowed 
by symmetry both above and below the structural phase transition at 395 K [9] whereas the existence 
of a finite intra-dimer DM interaction is allowed only below 395 K when the mirror symmetry of the 
(ab) plane is just slightly broken. Additional terms may be responsible for this disagreement, possible 
candidates being in-plane components of the inter-dimer DM interaction, which should also be allowed 
below the structural phase transition. 

We discuss now the issue related to the existence of a magnetic mode below the spin gap value [18, 21] . 
In order to reproduce the upward dispersion with fields B || c of the Tit, mode we had to choose a value 
of x which is greater than 0.5, otherwise this excitation would have displayed a downward dispersion. 
From Fig. 11 we observe that a ratio of x = J2/J1 > 0.5 implies that the position of the bound singlet 
state is below A. We suggest that this state is responsible for the observations of the 21.5 cm" 1 mode in 
Refs. [18, 21]. The presence of this excitation will also influence the specific heat measurements and, in 
conjunction with the finite intra-dimer interaction d a b, also the low temperature magnetization data whose 
quantitative understanding has not been achieved yet [1]. The existence of this magnetic mode is at odds 
with theoretical predictions [11, 12, 14]. However, the experimental finding of the set of A2 symmetric 
modes (T2 to T5 in Fig. 2, all of them below the two-magnon continuum starting at 2A w 48 cm -1 
and also not predicted by theory) already shows that the understanding of the spin dynamics in higher 
particle sectors is not complete. 
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Fig. 13. Energies (panels a and b) and intensities (panels c and d) of the spin excitations for B || c (right) and 
B _L c (left) from Fig. 3a-b. Symbols represent experimental points, solid lines are the results of 4-spin cluster 
diagonalization as described in the text, dashed lines are guides for the eye. 



Fig. 14 shows two low temperature Raman spectra taken in (a'b') polarization with two incoming 
laser frequencies, lo l = 1-92 and 2.6 eV. The point we make is that we observe two types of behaviors. 
Firstly we notice that the intensities of the group of A2 symmetric modes, Tie, and T2 is about the same 
for the two photon energies used. On the other hand, the modes corresponding to the group formed by 
the T e , Si, S2 and S3 resonances are more than two orders of magnitude stronger for lo l = 2.6 eV than 
for ujl — 1-92 eV. Other modes in the 50 to 70 cm -1 energy range also become visible in the uj l = 2.6 eV 
spectrum. 




Raman shift (cm 1 ) 



Fig. 14. T = 3 K Raman data in (a'b') polarization for ujl = 1-92 (filled squares) and 2.6 eV (empty circles) 
excitation energies. The data points corresponding to the resonantly enhanced singlet bound state at 29 cm -1 are 
multiplied by 0.15 and the corresponding line represents a Lorentzian fit. 
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Our data prove that the coupling to these two groups of excitations takes place by two distinct light 
scattering mechanisms. The resonance of the Bi symmetric T e magnetic mode, enhanced for lu^ — 2.6 eV, 
is similar to the one corresponding to Si, S 2 and S 3 excitations as well as to the behavior of the new 
modes seen around 55, 59 and 68 cm -1 . On the account of this similarity alone, one cannot identify 
the latter as magnetic bound states as well. However, the lack of energy shift in magnetic fields and the 
results of perturbational analysis regarding energy scales and symmetries [11] are not contradicting their 
interpretation as Raman active collective 5 = magnetic bound states [20]. 

We discuss below the nature of the two light coupling mechanisms to magnetic excitations. As for the 
set of A 2 symmetric modes we propose that the coupling takes place via the spin-orbit coupling which can 
be written in an effective form as (e in x e out )S t ^ t [22]. This interaction Hamiltonian probes excitations 
with A 2 symmetry and the calculated magnetic field dependence in Fig. 13 is also in agreement with 
the experimental results. The coupling to the Ti e and Tij modes from Fig. 10 can be understood if we 
invoke the usual effective spin interaction corresponding to the photon induced spin exchange process 
X^<i j>( e in' r ij){&out' r ij)Si'Sj- Here the sum runs over pairs of lattice sites, Si and Sj are the exchanged 
spins on sites i and j respectively, while is the vector connecting these sites [22]. Writing down the 
explicit expression of this interaction for several polarizations in the 4-spin cluster approximation (Figs. 9 
and 11) one indeed gets finite coupling in and B 2 channels for the triplet T le and T X j states. This 
explains the presence of the 25.6 cm" 1 modes in all polarizations except (RR). 

The difference in the coupling strengths seen in Fig. 14 is thus understandable because these two 
coupling mechanisms need not be simultaneously in resonance with the same high energy excited electronic 
states. The photon induced spin exchange Hamiltonian has been usually invoked in order to explain 
Raman active 5 = two-magnon type excitations in various magnetic systems [22]. SrCu2(B03) 2 is an 
example where this Hamiltonian, in the presence of singlet-triplet mixing DM interactions, can be used 
to account for coupling to 5 = 1 states. A remaining question is why do we not see the bound singlet 
mode with an energy below A for any of the two excitations used? In principle the photon induced spin 
exchange, resonant in this case at higher photon energies, could provide coupling to this excitation also, 
as to the 28.9 cm -1 mode denoted by Si. One possible explanation to be explored in more detail from a 
theoretical point of view is that the Raman form factors for exciting a pair of magnons both at k = is 
vanishing as opposed to the case of a pair of zone boundary modes. For instance this is the case when the 
Raman vertex is calculated for the 2D square lattice within the spin-wave approximation and using the 
Fleury-Loudon interaction Hamiltonian. Consequently, the 21.5 and the 28.9 cm -1 excitations could be 
both attributed to 5 — bound states originating from different parts of the reciprocal space and having 
substantially different binding energies. 

3.3 Summary 

We showed in this chapter low temperature Raman data on phononic and magnetic excitations in 
SrCu2(B03)2 • Regarding the former, in the to 350 cm -1 range we find several pairs of quasi-degenerate 
modes which have different symmetries. Group theoretical analysis suggests that the existence of these 
modes is related to quite different atomic vibrational pattern, i.e. in-plane and c-axis motions and as a 
result a quantitative investigation would be very interesting. Collective magnetic excitations were stud- 
ied in terms of symmetry, resonance and coupling mechanisms in zero and applied magnetic fields. The 
analysis of the 4-spin cluster shown in Figs. 9 and 11 allows us to understand the group symmetries of 
the zero field Brillouin zone center spin gap branches around 24 cm -1 confirming the picture of local 
elementary one-triplet modes. By considering an additional intra-dimer DM interaction we are also able 
to understand the observed selection rules and intensity variations of the spin gap branches in external 
magnetic fields applied parallel or perpendicular to the dimer planes. These selection rules also require 
that the energy of the 5 = two-triplet bound state made out of magnons confined within a unit cell 
is below A (in the 4-spin cluster this is equivalent to x > 0.5 in Fig. 10) suggesting a very high binding 
energy for this two particle excitation. The 4-spin cluster analysis fails to account for the two-triplet 
excitations which shows that they have contributions from the different parts of the reciprocal space. 
The existence of a set of four modes below the onset of two-triplet continuum, at 37.5, 40.8, 44.5 and 
50.9 cm -1 , in the A 2 symmetry channel, shows that further theoretical analysis is required in order to 
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understand the nature of these composite excitations. Finally, we identified two effective magnetic light 
scattering Hamiltonians responsible for the coupling to the magnetic modes which allowed us to explain 
their resonance behavior. 
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